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Abstract: We discuss rare kaon decays that violate charged- lepton flavour conser- 
vation in supersymmetric theories with and without R parity, in view of possible ex- 
periments using an intense proton source as envisaged for a neutrino factory. In the 
Minimal Supersymmetric Standard Model, such decays are generated by box diagrams 
involving charginos and neutralinos, but the limits from fx — * erf, fx-e conversion and 
AniK constrain the branching ratios to challengingly small values. However, this is 
no longer the case in i?-violating theories, where such decays may occur at tree level 
at rates close to the present experimental limits. Within this framework, we obtain 
bounds on products of LLE and LQD operators from the experimental upper lim- 
its on K° — * [r^e^ and K^' — * 7r =l= ' / u = ' : e =F decays. We also note the possibility of 
like-sign lepton decays — » ir^l^t^ 1 in the presence of non-zero 6l~^_r mixing. We 
conclude that rare kaon decays violating charged-lepton flavour conservation could be 
an interesting signature of R violation. 
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1 Introduction 



The recent Super-Kamiokande data have triggered an upsurge of interest in 
extensions of the Standard Model with massive neutrinos and/or violation of the 
charged-lepton numbers in processes such as \i — > ej, /i — > 3e, r — > ^7 and fi — > e 
conversion on heavy nuclei [@, |3], [|, || ||] . The present experimental upper bounds 
on the most interesting of these processes are: 



BR(n -> erf) < 1.2 x 10" 11 
BR(fi + -> e+e+e") < 1.0 x 10" 12 
i2(//~Tz -> eTi) < 6.1 x 10" 13 
BR(t -> ^7) < 1.1 x 1(T 6 



[7] 



Our main interest in this paper is in rare kaon decays, for which the current 
experimental bounds are: 
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Projects are being proposed that could be used to improve these limits signif- 
icantly, e.g., some of the powerful proton sources being proposed for neutrino 



factories |2T[ could provide intense secondary kaon beams. 

Any observable rate for one of these processes would constitute unambiguous evi- 
dence for new physics. The rates for such processes remain extremely suppressed 
if we simply extend the Standard Model to include right-handed neutrinos, but 
larger rates are possible in more ambitious extensions of the Standard Model. 
Supersymmetry is one example of new physics that could amplify rates for some 
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of the rare processes (p], ||D, either in the minimal supersymmetric extension of the 
Standard Model (MSSM) or in its modification to include violation of R parity. 

In previous works, rare charged-lepton decays and anomalous fi — > e conversions 
on heavy nuclei have received considerable attention ||, leading to a good under- 
standing of the correlations between the predicted rates and possible violations 
of leptonic universality and/or exotic Yukawa couplings with AL ^ 0. It is nat- 
ural to include rare kaon decays in this analysis, because strangeness-changing 
decays occur in the Standard Model. As mentioned above, there is now con- 
siderable discussion of intense proton sources to be used as muon sources for 
neutrino factories ||21|| . If the protons have high energy above ~ 15 GeV, as in 
some neutrino factory designs, they would also be copious sources of kaons. These 
might provide a new opportunity to study rare K decays with high statistics, and 
we are interested to know whether these might cast additional light on neutrino 
masses and mixing. Specifically, as we show in this paper, rare K decays that vi- 
olate charged-lepton number allow one, in the context of the MSSM, to correlate 
the (s)quark and (s)lepton mixing, whilst in i?-violating supersymmetry one can 
probe interesting products of Yukawa couplings. 

This paper is organized as follows. In Section 2, we analyze the decays K° — > 
yU ± e =F in the MSSM, finding rates that are relatively small, although for certain 
model parameters they may be within the reach of an imaginable kaon beam 
at a neutrino factory. In Section 3, we analyze the same decays in i?-violating 
supersymmetric models, finding that the rates could in principle be much larger, 
close to the present experimental limits. Indeed, these existing upper limits on 
these decays already impose interesting upper limits on products of LLE and 
LQD couplings. We also point out that b^-bji mixing in the presence of R 
violation could lead to the like-sign- lepton decays K ± — > n^^^, although the 
existing experimental limits on these processes do not impose interesting upper 
limits on couplings. Taken together, our results indicate that rare K decays could 
provide interesting signatures for supersymmetric models, in particular those with 
R violation. 



2 Rare Kaon Decays in the MSSM 

In this Section, we evaluate the rates for rare kaon decays in the MSSM with 



massive neutrinos, using the see-saw mechanism ]zz| , which we consider to be the 
most natural way to obtain neutrino masses in the sub-eV range. In particular, 
we assume Dirac neutrino masses of the same order as the charged-lepton 
and quark masses, and heavy Majorana masses M VR , leading to a light effective 
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neutrino mass matrix of the form: 



rn 



eff 
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(3) 



Neutrino-flavour mixing []23] may then occur through either the Dirac matrix m„ 
and/or the Majorana mass matrix M VR , which may also feed flavour violation 
through to the charged leptons. In non-supersymmetric models with massive 
neutrinos, the amplitudes for charged-lepton-flavour violation are proportional to 
inverse powers of the right-handed neutrino mass scale M VR , and thus the rates for 
rare decays are extremely suppressed ||. On the other hand, in super symmetric 
models one must also take the dynamics of sneutrinos v into account, and these 
processes may only be only suppressed by inverse powers of the supersymmetry- 
breaking scale, which characterizes and is at most 0(1) TeV ||. 

The magnitudes of the predicted rates depend on the details of the masses 
and mixings of sparticles including the sneutrinos. If their soft supersymmetry- 
breaking masses are non- universal at Mgut-, large rates are in general predicted |24| . 
However, even if the soft supersymmetry-breaking masses are universal at Mgut, 
renormalization effects in the MSSM with right-handed neutrinos spoil this di- 
agonal form 0, ||, [5] at lower scales. This is because the Dirac neutrino and 
charged-lepton Yukawa couplings and masses mi cannot, in general, be diagonal- 
ized simultaneously. Since both these sets of lepton Yukawa couplings appear in 
the renormalizat ion-group equations, the lepton and slepton mass matrices also 
may not be diagonalized simultaneously at low energies. 

To illustrate this point, let us consider the simplest example of a model based on 
Abelian flavour symmetries and symmetric mass matrices p5| |. Requiring large 
(2-3) mixing in this model [2B| severely constrains the possible flavour charges 
and thus the forms of the charged-lepton and the neutrino mass matrices. A 
representative example is given by Ansatz A of M: 
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where e = y/e = 0.2. We already see that the two matrices cannot be simultane- 
ously diagonalized and indeed, 
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Within this general framework, there is ambiguity in the specification of numerical 
coefficients in the matrix elements, which are expected to be of order unity. We 
return later to this point. 
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For a generic texture where the charged lepton and neutrino matrices are not 
simultaneously diagonal, the slepton mass matrix acquires non-diagonal contri- 
butions from renormalization at scales below Mqut- In the basis where mi is 
diagonal, these corrections take the form 0: 

^ K li^ 3 + a2 ' )ln ~^ :A -° ADm ^ /2 ' ^ 

where Ad is the Dirac neutrino Yukawa coupling, Mn is the scale where the effec- 
tive neutrino-mass operator is formed, a is related to the trilinear mass parameter 
Ai = am 3 /2, and m 3 / 2 is the presumed common value m of the scalar masses at 
the GUT scale. 

In the case of non-universal soft masses, these corrections are generically negli- 
gible. However, the rates for AL ^ processes are generally too large in such 
non- universal models ||. On the other hand, models with scalar-mass universality 
at Mgut, such as no-scale P7| and gauge-mediated models may yield accept- 



able rates for AL ^ processes. In such models, the contributions (||) related to 
neutrino masses dominate and lead to non-negligible rates for the lepton-flavour- 
violating processes which are determined by the off-diagonal terms in the Yukawa 
matrix A^. The various different solutions of the solar neutrino deficit |30j , 
with a small/large mixing angle and with eV or « 0.03 eV neutrinos, predict in 
general different rates for charged-lepton-flavour violation: the larger the \l — e 
mixing, and the larger the neutrino mass scales that are required, the larger the 
rates. Thus, degenerate neutrinos with bimaximal mixing may be expected to 
yield significantly larger effects than, for instance, hierarchical neutrinos with a 
small vacuum mixing angle. Note, in particular, that the just-so solutions to the 
solar neutrino problem with 5m 2 ~ 10~ 10 eV 2 predict, in the case of hierarchical 
neutrino masses, rates that are small, even if the first /second-generation neutrino 
mixing is large. 

In this MSSM framework, rare kaon decays are generated by box diagrams in- 
volving chargino and neutralino exchanges. For instance, for K° — > /i ± e =F we 
have the diagrams Q of Fig. 10. It is clear that K ± >° — > 7r ± ' /i ± e :f decays can be 
generated in a similar way, but in this case the experimental bounds are worse 
by almost two orders of magnitude, and we do not discuss them in detail. 

1 The contribution from the first diagram when the neutralino is a purely photino state has 
been discussed in 

2 Kaon decays in left-right symmetric models have been analysed in [j32| . In principle, using 
the mixing in the quark sector - in particular between s and d quarks - we can also generate 
K — ► fie by penguin diagrams. In the cases that all the quark mixing is in the down sector, or 
the right-handed mixing in the down sector is much bigger than the one in the left (which is 
the one bounded by Vqkm) the rates might be of relevance. However, we do not address here 
this model-dependent issue. 
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Figure 1: MSSM box diagrams for K° — > /j, e T . There is another neutralino 
exchange diagram corresponding to the permutation of the fi and e. Since x° is 
a Majorana spinor, there are contributions from the neutralinos that differ in the 
number of mass insertions. 



Our procedure for evaluating these contributions is as follows: 



• We first find the maximal squark mixing that is allowed by the neutral-kaon 
mass difference Arrix fl33| . 

• We next find the maximal slepton mixing allowed by \i — > e^j and fi-e 
conversion in nuclei in a model-independent way. 

• Having fixed these values, finally we calculate the rates for rare kaon decays. 



As we noted previously, the \x—e mixing is constrained by the form of the neutrino 
textures and thus by the recent neutrino data: in general, degenerate neutrinos 
with large angle MSW oscillations require smaller soft masses to be consistent 
with the observed rates. However, as we found in ||, even for the small-angle 
MSW solutions of the solar neutrino deficit, we can obtain large rates for values of 
m and well below 500 GeV. The latest Super-Kamiokande data []33J on solar 
neutrinos favour large mixing angles, which might suggest larger \x — e flavour 
violation. Thus, considering models with small mixing angles is conservative. 

For illustration, we focus here on one such model, namely variant A\ of the 
texture |5|), which has the numerical values of the ambiguous 0{1) coefficients 
fixed as described in |J. For definiteness, we choose its 'inverted' option with 
negative off-diagonal entries in the Dirac neutrino coupling matrix 0. Our results 
are presented in Figure 2, where we summarise our predictions for BR(K — > 
fie), in association with the predictions for BR(fi — > e^y) and BR(fi — e) for 
different values of tan/? and m^. We parameterize the supersymmetric masses 
in terms of the universal GUT-scale parameters m and m^, for sfermions and 
gauginos respectively, and use the renormalization-group equations of the MSSM 
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Figure 2: Illustrative predictions for BR(K — > fie), BR(fi — > ey) and BR(fi — e) 
for different values of 'tan (3 = 10 (left column), 20 (right column) andmi/ 2 = 150 
(dashed lines), 250 GeV (solid lines), as functions of m (in GeV). 
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to calculate the low-energy sparticle masses. Other relevant free parameters of the 
MSSM are the trilinear coupling A (for which we start with the initial condition 
A = —mi/2), the sign of the Higgs mixing parameter /1, and the value of tan/?. 
Models with different signs of \x give similar results: here we assume fi < 0. 

As expected, the larger the value of tan (3 and the smaller the soft supersymmetric 
terms, the larger the branching ratios, apart from certain cancellations. In the 
case tan/3 = 10 and 7774/2 = 250 GeV, we see that, for the range mo > 170 GeV 
where BR(fi — > ej) and BR(fi — e) conversion are consistent with the current 
experimental bounds (fj), BR(K — > fie) is at most 2 • 10~ 18 . However, for the 
same value of 777,1/2, when tan (3 = 20 we find a significantly larger branching ratio 
at small values of mo ~ 170 GeV. Moreover, for smaller mi/ 2 = 150 GeV, we gain 
almost two orders of magnitude when we consider mo in the low-mass window 
between 100 and 150 GeV. We recall that these lower values of m 1 / 2 ,m are 
consistent with accelerator constraints and generically yield cold dark matter 
densities in the range preferred by cosmology |35| . 

We do not discuss here other model textures for the mass matrices. Rather, our 
point here has been to demonstrate that, despite the limits from \i —>■ ej, fi-e 
conversion and Am^-, the branching ratio of K — > fie may be within the reach 
of the next generation of experiments, namely in the range 10~ 16 — ► 10~ 18 , at 
least if tan f3 is large and the soft supersymmetry-breaking terms are small. The 
sensitivities (§]) to fi — »■ ej and /x-e conversion could be improved significantly 
even before a neutrino factory comes into operation, and such a machine would 
offer enhanced prospects for probing them. It is therefore likely that the best 
prospects for discovering charged-lepton flavour violation may be offered by fi — e 
conversion and fi — e^. However, rare kaon decays provide a complementary probe 
which also gives information on the squark mixing, in the context of the MSSM. 



3 Rare Kaon Decays in R- Violating Supersym- 
metry 

We now discuss kaon decays violating charged-lepton flavour beyond the con- 
text of the MSSM. As is well known, the gauge symmetries of the MSSM allow 
additional dimension-four Yukawa couplings, of the form 

XLiLjEk, A LiQjDk, X"UiDjDk 

where the L(Q) are the left-handed lepton (quark) superfields, and the E,(D, U) 
are the corresponding right-handed fields. If all these couplings were present 
simultaneously in the low-energy Lagrangian, they would generate unacceptably 
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fast proton decay. Therefore, extra symmetries must be invoked to forbid all 
(imparity fl36|), or subsets (baryon and lepton parities |nj) of these couplings. In 
the latter case, very interesting baryon- and lepton-number-violating processes 
may occur 



Imposing electroweak SU(2) and colour SU(3) invariance, one finds that there 
are just 45 i?-violating couplings in total. Besides proton decay, there are many 
experimental constraints on these couplings, both individually and in various 
combinations, from the non-observation of modifications to Standard Model pro- 



cesses and of possible exotic processes |39], |40[ . In order to understand the possible 
hierarchies of /^-violating couplings, models of flavour symmetries have been in- 
voked For instance, it was found in previous work [53] that theories with 
symmetric fermion mass textures lead to the expectation that .R-violating cou- 
plings are small: < 10~ 3 -10 -4 for 100 GeV sfermion masses, whilst, in models 
with asymmetric fermion mass textures, dominance by a single coupling may be 
permitted, without however excluding several products of couplings from being 
non- negligible. 

In this paper, therefore, we allow the general possibility that several i?-violating 
operators may be large, and discuss the limits on their combinations that are 
obtainable from kaon decays. In this class of models, whilst fi — > ej occurs at 
the one-loop level, \x —>■ 3e, /i — e conversion, K° —>■ /i ± e =F and K ±,Q — > ir ±,0 fi ± e T 
may occur at tree level via different combinations of couplings. For instance, in 
the case of LQD couplings, \i — > 3e gives the limit |39| 



{LtQiD^LzQiDj) < 10" 4 
whilst /i-e conversion in Titanium gives: 

{L 2 QiDk){LiQiD k ) < 10- 8 
(L 2 Q j D 1 )(L 1 Q j D 1 ) < 10- 8 



in 



100 GeV 



(7) 



m 



f 



100 GeV 



m 



f 



100 GeV 



and the bound 



{L i Q 1 D 2 ){L i Q 2 D 1 ) < 10- 



m 



f 



100 GeV 



(9) 



(10) 



is obtainable from Amx 40 



What is the connection of these results with neutrino masses? If R parity is vio- 
lated, neutrino masses are generated via one loop diagrams involving the vertices 
VidkdjjVfdjdl for LQD operators and the vertices Vie k lj,vllje* k for LLE opera- 
tors. Focusing on LQD couplings and assuming that the left-right squark soft 
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mixing terms are diagonal in the physical basis and proportional to the associated 
quark mass, the induced masses are given by 



m Viil ~ " »J fc l '^ m , md I __z: — ^_ _| : — sl_ | (n; 



where f(x) = (xlnx — x + l)/(x — l) 2 , m^. is the down quark mass of the ith 
generation inside the loop, mj is an average of du and dm squark masses, and 
n c = 3 is the colour factor. 

Requiring that these contributions be consistent with the neutrino data gives 
bounds on the associated i?-violating products. We can see from the above 
expression that the heavier the fermions in the loop (including the associated 



fermion mass arising from the soft mixing term), the stricter the bounds |44 
For example, demanding m eA1 < 2.5 eV for sparticle masses of 100 GeV leads 
to A' 133 A233 < 3.8 • 10~ 7 , whilst for A 122 A' 2 22 the bound drops to 2.3 • 10~ 4 f44fl . 
For higher sfermion masses, larger /^-violating couplings are allowed. In Super- 
Kamiokande-friendly solutions with hierarchical neutrinos the bounds are stricter 
by two orders of magnitude. 

Note, however, that neutrino masses do not strictly constrain K — » fie (and in 
certain cases the rest of the flavour- violating-processes), since: 

• Neutrino masses may only constrain products of LLE or LQD operators, not 
mixed LLE-LQD products. 

• Even for the diagrams with products of only LQD operators, rare kaon decays 
involve quarks of the lightest and second-lightest generations. In this case the 
bounds from neutrino masses are significantly weaker, and the stricter limits 
come from the current measurements of the rare kaon decays themselves. The 
same is true for \i — > e conversion and even for \i — > via LLE couplings. For 
fermions of the first two generations, the bounds from the lepton-flavour- violating 
processes themselves tend to dominate. 

Two-body K° decays to muons and electrons proceed via the diagrams shown in 
Fig. |. 

At the quark level, the effective Lagrangian for such processes has the form f\ [40| : 



3 k mi 



KjkKn (s R d L ) (tjjtkR) + AifejA-2! (s L d R ) (ij R (k L ) 



%^(^d R )(j- L rik L ). (12) 



3 Note that, for energies of the order of the kaon mass, both s- and i-channel diagrams yield 
contact interactions. 
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Figure 3: Quark/ sfermion diagrams involving R-violating couplings that yield 
two-body K° — > £ ± £ :F decays. 



The two different contributions from s-channel and t-channel diagrams put limits 
on different couplings ||40|| , as we discuss below. 

We have derived Feynman rules for the relevant effective kaon, pion and lepton 
interactions. Based on these Feynman rules, we have recalculated the important 
kaon decay processes, and update the limits on the products of .R-violating cou- 
plings using the present experimental limits (|2|). The diagrams of Fig. |3] lead to 
the the following effective Lagrangian for K°£ + £~ interactions: 



C 



J h 



2ml 



K 



(13) 



where F. 



K° 



m 



K° 



/k/(w s + rrid), m s + ~ 0.15 GeV is the sum of the current 
masses of the s and d quarks, and fx = 0.1598 GeV is the kaon decay constant. 
The value of F K o is related to the pseudoscalar < 0\srfd\K° >= — F K o matrix 
element, and is obtained from fx by using the Dirac equations for quarks. All 
QCD corrections are included in this phenomenological approach. In the follow- 
ing, we assume that the .R-violating couplings are real and that only one of their 
products in flT2D is non-zero. 



We have implemented the Feynman rules in the CompHEP package j|5 
the effective Lagrangian (fL3|), and have obtained the following results: 



using 



K°- 



J k 



647r rn~.rriK° 



m 2 K0 i 



A(m K o,m e m Ek ) (14) 



2567T m~.ra.K-o 



+ m 4 



(ra| -ra? ) 2 \ 

— }A(m K o,m £ mt k ) 



m 



(15) 
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where A(a, b, c) = J [a 2 — (b + c) 2 ] [a 2 — (b — c) 2 ] is the triangle function. As 
there exist two similar contributions in the s channel, coming from terms with 
different couplings KjkKi2 an d ^ikjKzii as seen i n the first line in (0), we give 
only one of them in fll4"D. 

We obtain the following nominal numerical results for K° — > £ + £~ decay via v 
exchange, and the corresponding limits on the AA' products (for the numerical 
results for Y K o^+i- we have used the nominal values A = A' = 1 and mo{u) = 
100 GeV): 



/100 GeV\ 



2 



V°-* e + M - = 1-57 x 10~ 12 GeV, X i2 iX' il2 x ) < 6.2 x 10 



A , 12A ; 2lX (12!L^) 2 < 6 . 2xl o-» 

V mp J 

r>^ e - M+ = 1.57 x 10~ 12 GeV, A a2 A- 12 x ( 10 ° GeV ) < 6.2 x 10~ 9 



V TYl ~ 

^ x ( ™L2*L\' < 6.2 x 10- 

m,z 



T K o^ e+e - = 1.72 x 10~ 12 GeV, A m A' 12 x ( W ° GeV ) 2 < 1.0 x lO" 8 

\ 771 ~ J 

V mo J 

r>^+ M - = 1-42 x 10- 12 GeV, \ l22 X' il2 x ( 10 ° GeV ) 2 < 2.6 x lO" 7 

2 

-7 



A 422 A: 21 x f 10 ° GeV Y < 2.6 x 10- 7 (16) 

V mo ) 

For K° — > £ + £~ decay via up-squark exchange we have the following limits: 

r>^ e+ „- = 1.61 x 10~ 15 GeV, A' A' x f 100 GeV \ < L9 x 1Q -r 

V rriu J 

T K o_ e -„+ = 1.61 x 10" 15 GeV, A' ^A' 2 x ( W ° GeV ) < 1.9 x 10~ 7 

V rriu J 

T K o^ e+e - = 8.25 x 10- 20 GeV, A' X A' x (" 10 ° GeV ) < 4.7 x 10 -s 

" ma J 



T K o^- = 3.19 x 10- 15 GeV, A 2a A 2i2 x ( 10 ° GeV f < 5.4 x 10" 6 (17) 



We have used in our calculations the decay width r exp (i ; r°) = 1.273 x 10 17 GeV 
and the experimental limits on K° — > £ + £~ decay widths shown in (0). 
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We now discuss the diagrams for 3-body kaon decays to pions and two charged 
leptons, of which there are two qualitatively different kinds: 

• The kaon may decay into a pion of the same charge, in which case the 
leptons in the final state must have opposite signs: K — > tt £^£ and 
K° — > 7r°£ =F £ /± . The corresponding diagram for the first process is shown in 
Fig. i 

• The kaon may decay into a pion with the opposite charge, in which case 
the leptons in the final state must have the same signs: K ± — ► ti t £ ± £ >± . 
Representative diagrams for this process are shown in Fig. |5|. This process 
involves two heavy virtual particles, the W boson and a down squark. One 
should note that decay width of this process is directly proportional to the 
mixing between the left- and right-handed squark states, denoted by b l 
and 6r, respectively. If there is no mixing, this same-sign-lepton process 
mentioned vanishes. One can expect sizeable mixing only for squarks of the 
third generation (and especially in the high-tan (3 region), which is why we 
have used b L R in the diagram. 




v 

Figure 4: Diagrams involving R-violating couplings that yield the three-body 
leptonic decays K + — ► n + £~£ + . 



The possibilities for — > n £^£ r± or K° —>■ 7T £^£ decay may be further 
subdivided into two groups. 

• Diagrams involving only squarks, via which experimental upper limits bound 
products of LQD operators. 

• Diagrams involving also sleptons, which yield bounds on products of LQD 
and LLE operators. 

The diagrams of Fig. |] lead to the the following effective Lagrangian for K + n~£ + £~ 
interactions: 
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K J Til. ~ L > ' * ' - 1 



fc 2 fll; 
+ 



^T*%Xi2 (PK+P.T (tun At) K + (p K )7r-(-p n ) } (18) 



where M s ~ 0.49 GeV is the constituent mass of the s quark. In the general 
case, the vector matrix element of the K — > 7r transition is parametrized by two 
momentum-dependent form factors /+(£) and f-(t): 

<0\sYd\K + (p K )n-(- P7T ) >= 

< n + ( P7T )\sYd\K + (p K ) >= (p K +P,rU(t) + ( P K-PvYf-(t), 

where t = (px — Pn) 2 is the squared momentum transfered to the lepton pair. 
The experimental data on semileptonic kaon decay are adequately described by 
a linear approximation for / + : 

/+(*) = /+(o) [i + a+(*K) 

where A + ~ 0.03 and /_ ~ —0.31 ± 0.15. In the case of exact flavour SU(3) 
symmetry, the following relations hold: /+(0) = 1 and /-(0) = 0. Due to the 
Ademollo-Gatto theorem, the first form factor is renormalized only at second 
order in the S'?7(3)-violating interactions, and is therefore expected to be close to 
unity. The Ademollo-Gatto theorem is not applicable to /_, whose value is close 
to zero, as follows from experiments and the Callan-Treiman relation: 

/-«) = f -f- - f + {ml) ~ -0.15 . 

In 

For our estimations, we set /+(0) = 1 and neglect /_. To estimate the scalar 
form factor, we use the relativistic quark model, which gives [ : 

< 0\sd\K+(p K )7T-{- Pn ) >=< n + (p w )\sd\K + (p K ) >~ -2M 8 , 

where we keep only the leading term and drop the momentum dependence, and 
M s ~ 0.49 GeV as before. In the approximation of unbroken SU(3) symmetry, 
the corresponding form factors for the neutral kaons are smaller by a factor of 
s/2. 

We can now estimate the decay rates of the charged kaons: 

r K +^+ejetW = 256^3 Z A ~ m 3 ' J ^(t,m tj ,mi h )A(t,m K+ ,m^) 



(f-ml-ml)^ (19) 
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r K+ ^ + i 7 ^i) - 2{)lS7 ., mj m . 



i u i nb K+ 



( m K+-"V+) 2 

J A(t : m ij: me k )A(t : m K+: m n +) 



X 



(m 2 K+ + m 2 n+ + m\ + m\ - t 2 ) F 1 (t) - F 2 (t) 
if ~ ™ 2 . ~ m \)- 2ml+m 2 K+ 



x 



l + \ + {t 2 /ml) 



2 dt 2 



(20) 



where F x (t) = q max (t) + q min (t), F 2 (t) = 2(q 2 max (t) + q ma x{t)q m in(t) + q 2 min {t)) /3, 



and 



qmaxit) 



At 2 
1 

4t2 



(m 2 K+ -ml+ +m%. - m\^ - (A(t,m^,m 4 ) - A(t,m x +,m^+ 
\m 2 K+ - ml + + ml - m\$ - (A(f , m^, m 4 ) + A(t, m x+ , 



For the decay rates of the neutral kaons, the 577(3) factor of 1/2 and the corre- 
sponding masses for the neutral mesons m K o and m n o must be taken into account. 

We obtain the following nominal numerical results for K + — > n + l + £~ and K° — > Tr°£ + t 
decays via v exchange, and the corresponding limits on the AA' products (for the 
numerical results we have used the nominal values A = A' = 1 and mo{u) = 
100 GeV): 



1.38 x 10 n GeV, A i21 A^ 12 x f 100 GeV ) < L5 x 10 " 



A , 12A . 1X (1^X) 2 < 4 .5 X1 0- 
100 GeV\ 2 



r K o^ n o e+u - = 7.71 x 10~ 16 GeV, A m A' x ( —) < 7.2 x 10~ 6 

V mo J 



Aii2A- 21 x 



/100 GeVV 
V mo J 
100 GeVV 



< 7.2 x 10~ 6 



1.38 x HT 15 GeV, A a2 A' 12 x ( —) < 1.0 x 10 

V mo J 



Ai2iA- 21 x 



7.71 x 10~ 16 GeV, A a2 A^ 12 x 



100 GeVV 



mo 



< 1.0 x 10" 



100 GeVV 



m,~ 



< 7.2 x 10" 



, ,100 GeV\ 2 
A i2 iA^ 21 x ( ] < 7.2 x 10- 6 



^ 

mo J 
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2.14 x 1(T 15 GeV, A ?11 A' x f 100 GeV ) < 8-8 x 10 -5 



V 



, w ,'100 GeV\ 
AaiA' 21 x ) < 8.8 x 10" 5 



1.17 x 10~ 15 GeV, A 5ll A:, 9 x (" 10 ° GeV \ < 6 8 x 10 



K°^n°e+e- = X.X/ X XU beV, AniA a2 



/100 GeV\ 
V mo J 



2 



AiiiA- 21 x ( ) < 6.8 x 10 



-6 



= 7 - 58 x 10_16 GeV > *™Ki2 x f 10 ° GeV ) < 8-2 x 10 



, w /100 GeVV 

A i22 A 21 x < 8.2 x 10- 5 

V mo J 



Vo^mV- = 4 - 38 x 10_16 GeV - ^Kn x ( 10 ° GeV ) < 1.2 x 10~ 5 



/100 GeV\ 
\ rrir, J 



2 



A i22 A; 21 x ( ) < 1.2 x 10- 5 (21 



For the decays via up-squark exchange we have the following limits: 
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T K+ ^ +e+u - = 1.61 x 10- 16 GeV, A' 2ll A' H2 x f 1Q ° GeV Y < 1.3 x 10~ 5 
T K o^ 0e+u - = 9.06 x 10" 17 GeV, A 2il A' M2 x ( W ° GeV ) < 2.1 x 10~ 5 



r K +^ n+e - u + = 1.61 x 10~ 16 GeV, A'uiAU x ( 10 ° GeV ) < 3.0 x 10~ 6 

\ rriu J 

Y K ^ e - u+ = 9.06 x 10~ 17 GeV, A' X A' 2 x (^ 10 ° GeV ^) < 2 .i x 10~ 5 

V m u J 

V K+ ^ +e+e - = 2.38 x 10~ 16 GeV, A' H1 A' H2 x f 100 GeV ^ < 2 .7 x 10" 4 



r x o^o e+e - = 1.31 x 10~ 16 GeV, A' H1 A' li2 x [ 10 ° GeV ) < 2.0 x 10~ 5 

\ rn u / 



-ir ™ v w „ (™ CeVX 



2 



T K+ ^ +u+u - = 9.33 x KT 1 ' GeV, A' 2il A 2i2 x < 2.3 x 10" 4 

V m u J 

T K o^o u+u - = 5.48 x 10- 17 GeV, A' 2il A 2i2 x f 10 ° GeV ) < 3.4 x 10" 5 (22) 



We have used in our calculations the decay width T exp (K + ) = 5.314 x 10~ 17 GeV 
and the present limits on K + ix + l + £~ and K\ — > tt°£ + £~ decay widths @. 
Note that there are different constraints from the K + — > vr+e+jU^ and -fT + — > 
7r + e~yU + decays, because of the rather different experimental limits BR(K + — > 
7T + e + ri-) < 6.9 x 10~ 9 and BR(K + -> 7r+e-//+) < 2.8 x 10 _u . The limit obtained 



16 



from K° — > £ + £ is typically 1-2 orders of magnitude better than that derived 
from K + — > tt + £ + £~ decay. 




Figure 5: Diagrams involving R-violating couplings that yield three-body like-sign 
leptons decays K + — > ix^lfl^. 

As already mentioned, diagrams with non-zero sbottom-quark mixing may lead 
to like-sign leptons, as seen in Fig. |5]. They arise from the effective lepton-number 
violating contact interactions 

£ = - W ^ X jf LR (d-u iL ) ((J^fu qL ) + h.c, (23) 

" l d k 

where Vlr denotes left-right squark mixing matrix element. 

There are two different topologies, shown in the first and second row of Fig. |5|, 
respectively. Diagrams in the first row lead to an effective tree-level process, since 
the W^-boson and squark masses are much bigger then the typical energy scale, 
which is of the order of the kaon mass. Diagrams in the second row cannot be 
reduced to tree-level diagrams, since one has a neutrino propagator in the loop. 
However, this last set of diagrams give a contribution that is typically about 2 
orders of magnitude lower than the diagrams of the first row |46| |. Therefore, 



for the sake of simplicity, we neglect the second-row diagrams, and have derived 
an effective Lagrangian only for the first two diagrams of Fig. |5|. Two kinds of 
effective interaction appear: Standard-Model-like Kip and niv interactions and 
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new effective interactions related to i?-violating operators of the forms Kl v and 

The effective Lagrangian for those interactions take the following forms: 

= V us V2G F f KP » K {WIiAl) K + ( P k) (24) 
C n+i - Pt = V ud V2G F Up% (WZjAl) ?t + (k) (25) 

where fx and f w = 0.1307 GeV are the kaon and pion decay constants, respec- 
tively, G F = 1.1663910~ 5 GeV' 2 is the Fermi constant and V us ,V ud are CKM 
matrix elements, and: 

C K+{iJ) c Ui (d k ) = Kk2 *if LR F K+ (FFF^l) K+{ Vk ) (26) 

4 

t^yo^dk) = Kkl ^f LR F^ (I^F^l) kHp,) (27) 

4 

where F K + = m 2 K+ f K /(m s + m u ), F w + = ml+f n /(m d + m u ), m s + m u ~ 0.15 
GeV, and m d + m u ~ 0.01 GeV. 

The Feynman diagrams for K + — > ir~£f£j decay are given in terms of these 
effective interactions, as shown in Fig. || They yield the matrix elements: 

M = KklX '^i RF * + V^f*™* ((^F P f^ K + { Pk )k-^) (28) 

and 

M = ^ 2X ^l RFK+ V^Urn^ {W^f P fi iR ) K+fafr-M. (29) 

Recalling that V us /V ud ~ 0.2, F K + ps F w + and fx ~ f n , the main contribution 
comes from the last matrix element with i — 2, because the matrix element is 
chirally suppressed in the other case. As an example, one can obtain a constraint 
on the product of A / 2A . 2 A / llfc and Vlr from K + — ► 7r~e + / u + decay. In this case, we 
have the following numerical result []: 

T(K+ - ,TeV) = Vl R {\' 2k2 \' llk f x [lOOGeVV x 2 . 5 x 10 -« GeV (30) 

V m 4 / 

or 

^(A' 2fc2 A' llfc ) x (±°1^) 2 <10. (31) 
V m 4 / 

4 In a similar constraint was found with a different choice of diagrams. 
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It is apparent that kaon decay into a pion and a like-sign lepton pair is too strongly 
suppressed to be useful at present: the corresponding bounds on the A' A' product 
are currently of the order of 100. 



Figure 6: Diagrams for the like-sign lepton decay K + — » n~£f£j, in terms of the 
effective Standard- Model-like interactions Ktv and ix£v and effective Kt c v and 
nl c v interactions related to R-violating operators. 

4 Conclusions 

We have discussed in this paper flavour-violating decays of kaons into charged- 
lepton pairs in supersymmetric theories, in both the Minimal Supersymmetric 
Standard Model and .R-violating models. In the first case, these decays are 
generated by box diagrams involving charginos and neutralinos, and both the 
squark and the slepton mixings enter in the analysis. The process looks promis- 
ing for correlating the (s)-quark and (s)-lepton mixing by a combined study of 
rare charged lepton and kaon decays. Despite the limits from fi — > ery, /i-e con- 
version and Arrix, the kaon decay branching ratios for large tan/? and small 
soft supersymmetry-breaking terms may be accessible to a future generation of 
experiments using new intense proton sources. 

In the case of .R-violating supersymmetry, such rare kaon decays may occur at 
tree level. In this case, /i — > ej again occurs via one-loop diagrams, whilst /i-e 
conversion may also occur at tree-level, but via a set of operators different from 
those relevant to kaon decays. In this framework, we studied the expected rates 
for the decays K — > /i ± e =F and K —>■ rrfie, for all two- and three-body processes. 
Using the current experimental data, we obtained the bounds on products of 
LLE and LQD operators summarized in ([16]) and (|T^ ). We have also noted the 
possibility of like-sign lepton events in the presence of non-zero -bn mixing, but 
for this to occur at significant rate one would need large .R-violating couplings. 

Our final conclusion is that lepton-flavour-violating rare kaon decays have the 
potential to provide important information on the issue of flavour physics. Any 
future observation would, in addition, help distinguish between different super- 
symmetric theories. 



K+ 
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Appendix: The rate for K — ► fie decay via box diagrams in 

the MSSM 



BR{K - fxe) = ^ - 2n (\K L \ 2 + \K R \ 2 ) (32) 



The branching ratio for K — > fie is given by: 

2.65A 2 
2G F sin 2 9 W 

where A = (m^ + m s )m^/m 2 K , and K L , K R are given by the following expressions: 

K L = K c L + (K n L {1) + K n L {2) ) (33) 
K R = K% + (K^ + K^) (34) 

with 

T^c 1 t ( ^ n R{d) r ,R{d)* r ,R{t) ^Bid) n L{d)* n L{t) n R{£)*\ / Q r\ 

-^L — ^ J 4(A,B,X,y) I — — <-'dAJs: L/ sBX ^^By^eAy ' ^dAX^ sBX ^ fiBY^eAY I l J0 J 

1, ( r L{d) r R{d)* r L{l) r R(£)* , r L(d) r L(rf) V fiW r K(£)*\ 

K R = -KRu+r (36) 



where m^ A B and mr denote chargino and sneutrino masses, in the chargino 



b XA,B " b lx : 

contribution. Moreover 



, d k k . . 

iMa,b,x,y) = I ( 27r )4 ( fc2 _ M|J(fc 2 - M|J(A; 2 - m? )(k 2 - m? ^ (3?) 



74/ 



i d k X . . 

^4(A,s,x,y) = / ( 2 7T) 4 (A; 2 -M| A )(fc 2 -M| fl )(fc 2 -m 2 )(jfc 2 -m? v (38) 



(y - 



whilst the mixing coefficients C^' 1, appear in the fermion-sfermion-chargino in- 
teraction Lagrangian: 

+di(Cf^xPR + C\axPl)Xa^x 

+u 1 {C^ U x ) Pr + Cfj$P L )xtd x + h.c. (39) 



and their explicit expressions are given, for instance, in 
The fermion-sfermion-neutralino interaction Lagrangian is similarly written as 

Ant = fi(N$$P R + N^P L )x° A fx (40) 
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where / stands for I, u, d and u. The neutralino box contributions corresponding 
to the permutations of the fi and e in the external lines, K^(2) and K^(2), are 
then found to be 

K n(l) _ 1 j /A R( d ) N R(d)* N R(£) M R(£)* N R(d) *jL(d)* A jL(t) A jR(e)*\ /,-, \ 

K L - ^4(i,B,X,y) {-^^dAX^sBX ^fiBY^eAY + ^ dAX iV sBX iV fiBY 1 ^ eAY I l 41 J 

-^h(A,B,X,Y) m x° A m A {^dAX^sBX \iBY eAY + AI ^dAX^sBX \iBY eAY ) 

K n(l) = _ K f)\^ R (42) 

and 

K n(2) _ 1, ( ^ N R{d) N R{d)* N L{£) N L(i)* N R(d) A jL(d)* A jL{£) N Wfi*\ (A o\ 

K L - -^ J 4(A,B,X,Y) IT} d AX sBX ^ BY eAY dAX sBX fiBY eAY I l 4 ^J 



, 1 r „ „ ( AfL(d) Af R(d)* N L{i) A jR(i)* . Ar L(d) A jL(d)* AJ L{1) Af L(e)* 

+ ^h(A,B,X,Y) m x A m x° B {^dAX^sBX ^^BY^eAY + A^dAX^sBX ^fiBY^eAY 

K (2) = -Kf^ R . (44) 
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